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In this paper we apply second-order gauge-invariant perturbation theory to investigate the pos- 
sibihty that the non-hnear couphng between gravitational waves (GW) and a large scale inhomo- 
geneous magnetic field acts as an amplification mechanism in an 'almost' Friedmann-Lemaitre- 
Robertson- Walker (FLRW) Universe. The spatial inhomogeneities in the magnetic field are consis- 
tently implemented using the magnetohydrodynamic (MHD) approximation, which yields an addi- 
tional source term due to the interaction of the magnetic field with velocity perturbations in the 
plasma. Comparing the solutions with the corresponding results in our previous work indicates that, 
on super-horizon scales, the interaction with the spatially inhomogeneous field in the dust regime 
, induces the same boost as the case of a homogeneous field, at least in the ideal MHD approximation. 

This is attributed to the observation that the MHD induced part of the generated field effectively 
i-C ' only contributes on scales where the coherence length of the initial field is less than the Hubble 

scale. At sub-horizon scales, the GW induced magnetic field is completely negligible in relation to 
the MHD induced field. Moreover, there is no amplification found in the long- wavelength limit. 
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O ■ I- INTRODUCTION 

(N 

' Astrophysical observations indicate that almost all environments in the Universe are magnetised, and the more 
' we search for extragalactic fields, the more pervading they are revealed to be. Such cosmological magnetic fields are 
found in galaxy clusters, disk and spiral galaxies as well as in high-redshift condensations (see ^| for reviews). A 
fascinating and as yet unsolved question is how did these fields originate? The current properties of magnetic fields 
should, in principle, reflect their past and give clues to their origins, so we rely on observations to contribute to finding 
an answer to this important question. 

Faraday rotation and Zeeman splitting measurements indicate that galactic magnetic fields at high redshifts exist 
with roughly the same strength, 10"'' to 10~^ C, as those found in the Milky Way [^|^. The common properties of 
large scale fields in different galaxies indicates that their origins may be intrinsically connected to the cosmological 
repercussions of the interplay between gravitational and gauge interactions 0. This suggests that their origin may 
d ' be primordial, in which case their presence could be related to the physics of the very early Universe. One example 
is Big Bang Nucleosynthesis, where the interaction of a magnetic field with the magnetic moment of a neutrino may 
have given rise to a spin-flip and change of its handedness, introducing an additional neutrino degree of freedom 

Depending on their spectrum, magnetic fields existing in proto-galactic clouds with strengths of 10"^^ to 10~^ C may 
have played a significant role in structure formation [6j . The Lorentz force that acts on charges in an inhomogeneous 
(i.e. curls ^ 0) magnetic field has been shown to induce peculiar velocities 7] which seed density perturbations, 
and in so doing alter the gravitational instability picture. Moreover, hyper-magnetic fields, hypothesised to emerge 
during the electroweak phase, have also been identified as a possible source of the observed baryon asymmetry of 
the Universe 0. For these reasons, the determination of the origin and properties of cosmic magnetic fields is of 
extreme importance in cosmology. This makes magnetogenesis, the determination of a self-consistent theory for the 
generation of cosmological magnetic fields with the strengths and on the scales measured today, one of the 'hot' topics 
in modern cosmology. The most popular theories include the amplification of a small field by the galactic dynamo and 
the adiabatic proto-galactic collapse at the start of structure formation. Although these mechanisms are shown to 
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yield substantial enhancement, they are not self-sufRcient as they presuppose the existence of seed fields. In addition, 
these seed fields must satisfy very stringent strength and size criteria in order for the generated fields to agree with 
the magnitudes observed today. It follows that the problem we face is to provide a mechanism that induces a large 
enough amplification of a weak pre-existing seed field, so that the aforementioned mechanisms are physically viable. 

In this paper, we extend the work of Betschart et al. which investigated the coupling between a large-scale 
homogeneous magnetic field and the gravitational wave spectrum which accompanies most inflationary scenarios. 
This work built on earlier work by Tsagas et al. T^l in which the same interaction was studied within the weak field 
approximation 11, 12J. The analysis in ^ demonstrated that this coupling can lead to an amplification, provided 
the dimensionless shear anisotropy a/H at the end of inflation is larger than 10"**°. 

In this investigation, we consider the general case where the original magnetic field is inhomogeneous over a typically 
observed coherence scale. By comparing this analysis with our treatment of homogeneous fields in Q, we aim to 
determine the implications of placing restrictions (such as homogeneity) on the properties of primordial magnetic 
seed fields. 

The highly non-linear nature of the Einstein's field equations makes finding exact solutions as well as applying 
numerical techniques extremely difficult. In order to solve them analytically, severe symmetry assumptions are often 
required to simplify the physical models, which then restricts their applicability. For cosmological applications, using 
a perturbative approach, which entails decomposing the real physical Universe into a family of spacetimes, yields 
surprisingly good results. This allows us to encode the inhomogeneities we see today as perturbations expanded 
around a fictitious idealised background model, most commonly taken to be the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetimes. Before we can implement perturbation theory in a self-consistent way, we must give due 
attention to the issue of gauge- invariance which has historically plagued such studies |l3 • 

If we consider a weak large-scale magnetic field residing in a background FLRW model as done by Tsagas et al. |^ , 
the interaction with linearised gravitational waves manifests itself as a first-order perturbation of this background. 
The problem with this approach is that it is not gauge-invariant in a strict mathematical sense, due to the fact that 
the magnetic field introduces a preferred direction and therefore breaks the isotropy of the background FLRW model. 
This problem is partially overcome by assuming that the magnetic field is weak and that its contribution to the 
energy-momentum tensor is such that it does not disturb the isotropy of the FLRW background 12J. 

A completely self-consistent solution to this problem is obtained by treating any seed magnetic fields as a first order 
perturbation and including the interaction with gravitational perturbations by going to second order in perturbation 
theory JJ. 

The introduction of magnetic spatial gradients at linear order requires a more subtle treatment of the associated 
spatial currents, requiring the use of the magnetohydrodynamic approximation (MHD) to provide a framework in order 
to obtain a tractable solution. This single-component fluid model provides an accurate description of a two-species 
plasma when effects occurring over much larger time and length scales than those characteristic of plasma effects are 
studied. This allows one to handle low-frequency phenomena in a magnetised plasma using the standard machinery of 
fluid dynamics. This reduced description is achieved by deflning appropriate one-fluid variables representing the bulk 
quantities, while Ohm's law provides a consistent treatment of the associated electric field. It is common practice to 
take advantage of the high conductivity, cr, of the young cosmic plasma and employ the ideal MHD limit (ct oo). 
In this limit, the flux lines are effectively glued to the plasma fluid elements. In the non-ideal MHD limit [lj|, where 
a is assumed to be large but not infinite, the electric field enters at linear order in the case of a first-order magnetic 
field. Here we focus on the ideal MHD case which is shown to be equivalent to assuming that the observed electric 
field vanishes in the rest frame of the fluid. 

The mathematical framework we use is the 1-1-3 covariant approach [Tslll^lr^ to perturbation theory which allows 
Maxwell's and Einstein's equations to be written in an intuitive and simple fashion p^ . The covariant deflnition 
of the variables ensures that their connection to physically and geometrically significant quantities is immediately 
transparent and their exact presentation gives them meaning in any spacetime. Most importantly, the identification 
of gauge-invariant (GI) perturbation variables at a given order is relatively straightforward. 

It is found that the generated magnetic field consists of contributions from two sources. The first stems from the 
interaction between GWs and the 'background' magnetic field, the second comes from the rotation of the induced 
electric field, which is caused in MHD by the 'background' magnetic field interacting with the velocity perturbations 
in the plasma. Comparing the solutions with the corresponding results in Ji] reveals that, at super- horizon scales, 
the interaction of GWs with a spatially inhomogeneous magnetic field in the dust regime yields an amplification of 
the same order of magnitude as found in the case of a homogeneous field, at least in the ideal MHD approximation. 
The MHD induced part of the magnetic field becomes important only at sub-horizon scales, where the GW induced 
contribution is negligible. In other words, the contribution stemming from the GWs dominates over plasma effects in 
the ideal MHD limit at super-horizon scales, whereas the roles are interchanged at sub-horizon scales. 

The units employed in this paper are c = h = 1 and n = SnG = 1. 
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II. PERTURBATION SCHEME 

We employ the same perturbative scheme as in Betschart et al. Jj . Since the presence of a magnetic vector field 
in the FLRW background does not yield a strictly gauge-invariant system, we introduce such fields as a perturbation 
of the FLRW background. Similarly, the gravitational and velocity perturbations are required to vanish in this 
background for them to be gauge-invariant. Since we are interested in the interaction of magnetic fields with linear 
gravitational wave distortions, we need to treat this problem at second-order in perturbation theory. Using the 1-1-3 
covariant approach lid Il7| | , we expand the physically relevant variables in terms of two smallness parameters 
to distinguish between the magnitudes of the inhomogeneous magnetic field (~ e^) and the amplitude of the GWs 
(~ Cg). It follows that the magnitude of the interaction of interest is of order 0{eB^g)- Since we are only interested in 
the cross interaction of the magnetic and gravitationally sourced perturbations, we only retain mixed terms of order 
0{eBig) and neglect those of order 0{e^) and O(e^). In fact, such terms always appear in the calculations concerning 
the induced magnetic field, multiplied by terms of order ©(es) and therefore do not lead to inconsistencies in the 
perturbation scheme. 

The perturbation spacetimes may be split as follows p]: 

9 B ~ Exact FLRW as the background spacetime, ©(e*^); 

• J-i = Exact FLRW perturbed by an inhomogeneous magnetic field whose energy density and anisotropic stress 
are neglected, ©(es); 

• = Exact FLRW with gravitational and velocity perturbations [Mj 0{eg); 

• S — J-i + J-2 allows for inclusion of interactions terms of order 0{tB^g)- 

We will generally refer to terms of order ©(es) and 0{eg) appearing in T as 'first-order' and to the cross terms 
0{tB^g) appearing in S as 'second-order'. 

As we will discuss in more detail in the next section, the presented hierarchy of spacetimes given above is only 
justified if the electric field is at least of second order. However, we will assume the ideal MHD case, which indeed 
implies that the electric field must be second order. 

Before we turn to the MHD approximation, we first describe the basic equations describing the background and 
first-order spacetimes keeping the equations as general as possible. 

A. FLRW background 

In order to perform an 1-1-3 decomposition of any spacetime, we need to introduce a universal reference 4- velocity 
field u° relative to which all motion is defined and quantified. In accordance with the observed average recession 
of the galaxies, we assume that the matter in the Universe has a locally well-defined preferred motion that can be 
represented by a unique 4-velocity vector field satisfying Uau" = —1. Based on the Copernican principle, we can 
assume that this holds at each point in the Universe. We then introduce a family of observers, called the fundamental 
observers, travelling such that this field represents the congruence of their worldlines. In so doing, any observations 
made are those relative to this preferred frame. In choosing this field to coincide with the average velocity of the 
matter in the Universe, it acquires an 'invariant significance' such that the covariant quantities at every point, which 
are defined with respect to can be decomposed uniquely 0|. 

FLRW models are spacetimes that are spatially isotropic and homogeneous about every point. Relative to the 
congruence of fundamental observers with 4-velocity u", the kinematics are assumed to be locally isotropic. This 
implies that all tensorial quantities such as the acceleration vector iia = m'' Vf,Uo, the shear aab = D<aWf,> and the 
vorticity LOab = D[a ^6] must vanish to eliminate preferred direction in the spatial sections: 

= lia = O-ah = LOab, (1) 

The constraints 

TTab = (7a = , (2) 

ensure that the energy-momentum tensor and thus the Ricci tensor are isotropic and indicate that a perfect fiuid is a 
necessary requirement of this model. These restrictions mean that the electric and magnetic components of the Weyl 
tensor and hence the Weyl tensor itself, vanish identically 



Eab — Hab = =4' Cabcd — 



(3) 
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We can then infer that these models are conformally flat. Furthermore the spatial uniformity forces the spatial 
gradients of the energy density fi, the pressure p and the expansion Q to vanish 

= D,// = D,e = D,p. (4) 

As usual, the spatial derivative = hj^ V;, is obtained by projection of the spacetime covariant derivative onto the 
3-space (with metric hat = gab + UaUj,) orthogonal to the observer's worldline. As a consequence, the key background 
equations are the energy conservation equation 

A + e(/i + p) = o, (5) 

the Raychaudhuri equation 

e = -ie2-i(^ + 3p)+A, (6) 

and the Friedmann equation 

. + A4e' + f, (7) 

describing the intrinsic curvature of the homogeneous and isotropic 3-spaces. The curvature constant K indicates the 
geometry of the Universe and can be normalised to K = — 1,+1,0 for spatially open, closed and flat Universes. 



B. First-order perturbations 

1. The inhomogeneous magnetic field Ba 

We assume that the seed magnetic field Ba residing in the Ti spacetime is inhomogeneous over its typical coherence 
scale. The spatial gradients Di,Ba are thus of order ©(es). Given that the magnetic field is a first-order perturbation 
on the background, the magnetic anisotropy nab — —B<^aBb> ^ C>{e%) can be neglected ■ Since the associated 
electric field is perturbatively smaller than the magnetic field and enters only in S as will be argued in the next 
section, the magnetic induction equation has the form 

B<a> + iQBa = 0. (8) 

Thus the magnetic field decays as 

where a denotes the scale factor, i.e., 8 = 3 ct/a = iH and H denotes the inverse Hubble length. The adiabatic decay 
evident in equation ^ arises from the expansion of the Universe which conformally dilutes the field lines due to flux 
conservation . By taking the spatial gradient of equation Q , it is easy to show that the gradient of the magnetic 
field evolves as Di,Ba ^ a^^. Note also that the induction equation does not discriminate between homogeneous and 
inhomogeneous magnetic fields. 



2. Gravitational waves 



In the covariant approach to cosmology, linearised gravitational waves are purely tensorial and are monitored via 
the electric (Eab) and magnetic {Hab) Weyl constituents, which are not sourced by rotational (vector) and density 
(scalar) perturbations 19, 20]. The transverse (divergence-free) nature of these projected, symmetric, trace-free 
(PSTF) tensors means that we only need to eliminate their vector parts in order for them to characterise frame- 
invariant GWs. We isolate the linear tensorial modes by imposing the constraints 

= D,/i = Dap = D,e = u;a^ iia. (10) 

These restrictions ensure that the sources of vector modes (spatial gradients and vector perturbations themselves) 
vanish and lead to the constraints Hsl 



= DVafc = J^^Eab = D"i/afc = Hab - CUrl CTab- 



(11) 
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Since the shear tensor is coupled to Hab and Eab, it can also be used as a measure of gravitational waves. The 
tensorial gravitational waves are governed completely by a closed wave equation for the shear. At linear order, this 
wave equation is given by 

^<afc> - D^a^f, + f ecr<„f,> + (^e^ + - |p + |A) aab = 0. (12) 

Observe that in general the RHS of equation (|12(l is nonzero (see for the case of irrotational dust spacetimes) . 

III. BASICS OF IDEAL MHD 

Since the mean free path between the electron-ion collisions in a typical plasma can be short compared to the 
characteristic MHD length scale, it is not always obvious that a fluid description is indeed valid. On small scales 
these interactions are frequent and cause the two species to move relative to each other, generating charge separation 
effects referred to as plasma oscillations. If we consider much larger scales on which the individual collisions are not 
explicitly seen (i.e. low-frequency phenomena), the different species are observed to move together with a common 
average velocity, allowing an effective single fluid description of this two-component system. If, in addition, the 
characteristic MHD length scale is much greater than the plasma Debye length and the gyro radius, then MHD gives 
an accurate description of low-frequency phenomena in a magnetised plasma (see, for example, Refs. [T^ l22l|'). 

We now pay attention to how the currents, which are established due to the net motion of the fluid by induction, 
modify the field and in so doing couple the hydrodynamical equations to Maxwell's equations via Ohm's law. 

We perform the calculations in an irrotational Universe with vanishing cosmological constant A, and assume p = 
from now on since it has been argued that the reduced MHD description of a plasma is only valid in the cold limit 
(i.e. p — and non-relativistic motion of the plasma particles). These assumptions simplify the analysis tremendously 
while still allowing us to make a meaningful comparison with the results obtained in the homogeneous case for dust 

As a further simplification, we adopt the geodesic frame, in which the acceleration of the fluid frame vanishes to 
all orders. 



A. Covariant Theory 

We now turn to the evolution equations of magnetohydrodynamic variables. We assume charge quasi- neutrality of 
the plasma (i.e. the number densities of the electrons, rig, and ions, n^, are roughly equal such that the total charge 
Pc effectively vanishes, pc = —e{ne — rii) « 0). In the spirit of the MHD approximation, we choose to formulate the 
equations using the magnetic field as our primary variable and use Maxwell's equations to express the electric field 
and currents in terms of _B°. In order to identify the perturbative order of the terms which we can eliminate without 
excluding physically important effects, the exact equations that are necessary to fully describe this system are stated. 

1. Maxwell's equations 

Under the assumption of charge neutrality and vanishing vorticity, the behaviour of electromagnetic fields and 
currents in curved spacetimes is governed by Maxwells equations in the following form (see, e.g., [14. .16, .29]): 



B<a> + ^eBa = aabB'' - curl Ea, (13) 

E<a> + leEa - curl Ba-j<a>+<yabE\ (14) 

B-Ba - 0, (15) 

D°Sa = 0, (16) 



where the last two constraint equations are supposed to hold at all orders. 

2. Equations of motion 

We assume that the interactions of the ions and electrons collectively isotropise their motions, such that in a chosen 
frame the properties of the fluid on macroscopic scales can be described in terms of an average velocity w". We take 
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this mean motion to be the center of mass velocity of the electron-ion system, defined by 

Me + Mi 



(17) 



where u° coincides with the velocity of the fundamental observer at zeroth order, that is, v"' = relative to the 
fundamental observer at lowest order. At higher orders, the behaviour of the electromagnetic fields also depends on 
this bulk velocity. 

Ohm's law is generally formulated in the local rest frame of the conducting fluid and is assumed to hold at all 
orders: 

j<a> = <j{Ea + eabcV^B') , (18) 

where the conductivity a is taken to be a constant for simplicity. The 3-vector represents the field as observed 
from the rest-space of the fluid and the second term in (|18|l is the apparent electric field associated with the bulk 
velocity We make the standard assumption that the cosmic medium is infinitely conducting, and consequently 
apply the ideal MHD limit. This is valid considering the early epoch in which the interaction takes place and is also 
consistent with the treatment of a homogeneous magnetic field in a spatially flat Universe in Betschart et al. j^; 
the gravito-magnetic interaction was shown in section (IV. B) of |^ to generate the same magnetic field irrespective 
of conductivity. By letting cr ^ cxd in Ohm's law H18|l . we find that Ea + eabcV^B'^ must tend towards zero in order 
for the spatial current j^"^ = PeVe + Pi"^" to remain finite. The electric field is now determined jointly by the fluid 
velocity and the magnetic field: 

Ea = -eabcv''B' . (19) 

It follows that the electric field is at least of second order and vanishes at lower orders. At first order, the magnetic 
field B'^ determines the current via curl = j^"^, while at second order the current follows from equations (|14|l and 
(|19|l . Using charge neutrality, the evolution of the total energy density £ = + l^-i and the center of mass velocity 
follow from the corresponding total energy and total momentum conservation equations [l4l |29|| : 

£ + e£ = -VlaiSv"), (20) 

£ {v<a> + lOVa) = -£ {v'''DbVa + a^bf ') - EabcfB' , (21) 

which hold in the cold plasma limit in the geodesic frame. 



IV. SECOND-ORDER PERTURBATIONS: THE INTERACTION 



We look to Maxwell's equation to determine the nature of the interaction between GWs and the first-order magnetic 
field . If the back-reaction of the induced field with the shear and the center of mass velocity is ignored, the induction 
equation takes the form 

S<a> + feB, = (7„fcB'' + 2D^(t;[„Sfc]) , (22) 

where the second term on the RHS —curl Ea describes the dragging of the field lines by the fluid. We are faced with 
the problem of removing the primary magnetic field component of B" from the LHS of 122|) to ensure that it is truly 
second-order. Given that the magnetic spatial gradients are now retained in J-'i, the commutation relation used as an 
example in the homogeneous case in Betschart et al. 1^, is now consistently satisfied when the power series expansion 
ofS", 

B"^ eeB'l + egeeB^ + Oielel) , 

is applied. The reason why an inconsistency arose previously seems to stem from the requirement that the first-order 
magnetic field is homogeneous, which messed up the commutation relations. A similar inconsistency was encountered 
when studying the gravito-magnetic interaction in the vicinity of a Schwarzschild black hole, where staticity was 
imposed upon the first-order magnetic field (see [2jj| for details). It is interesting to remark, however, that constraining 
the magnetic field to be solenoidal does not pose any problems whatsoever. 

Although the above expansion does not immediately appear to be invalid, in the event that an inconsistency might 
exist, we choose to represent the magnetic field using the second-order gauge invariant (SOGI) variable 



Pa = B^a> + f Q-S, 



(23) 
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first identified in Betschart et al. We select the interaction variable la = cTabB'' and define the variable Fa = 
—curl Ea = 2D^ (^V[aBi,]j . We note that the electric field is of the order 0{eB£g) and thus enters the S spacetime. We 
can now restate Maxwell's equations in 5 as a system of differential equations in terms of these SOGI variables 

Pa = la+Fa, (24) 

B'^Ea = 0, (25) 
D"Ba = 0. (26) 

To close the system, we use the velocity propagation equation (|21|l . Given that both v'^ and B'^ are individually 
regarded as first-order, only the linear part of this equation is needed: 

V<a> + ^QVa - . (27) 

From equation (|24|l . we see that the generated magnetic field i?° can be found directly by integrating a linear 
combination of the /" and F° solutions once they are found. 

It is convenient to re-scale the magnetic field variable by defining Ba = Ba (a/oo)^. The time dependence of 
found in the final solutions then describes the evolution of the generated field relative to the 'background' field. The 
main variables become 

Pa^[^yB<a>, la^TabB", = 2D' (v^aB^]) . (28) 

Using H = Q/3 = a/awe can restate the important equation H24|) in terms of these variables as 

B<a> =Ia+J'a. (29) 

Integrating this equation with respect to proper time yields B"^; the constant of integration is determined by the 
physical condition that at the time Iq, when the interaction begins, there is no generated magnetic field, so we have 
B^^ Bg = initially 



V. EVOLUTION EQUATIONS FOR THE MAIN VARIABLES 

It is evident from equation (|29l) that the SOGI magnetic field can be extracted directly by integrating a linear 
combination of the solutions for the interaction X° and the electric field rotation !F°'. We turn to find the evolution 
equations for these variables, maintaining generality. 



A. Harmonics 

We employ the standard harmonic decomposition |l3l l25j | to deal with the Laplacian operator present in the wave 
equation of the shear. It is standard procedure to assume that the time and spatial dependence of each variable 
is separable, so that the variable can be expressed as the product of the time and spatial parts. This operation 
effectively decomposes the differential equation for the time evolution of a perturbation variable into separate equations 
describing the time evolution of each harmonic component which characterised by a comoving wavenumber k. Since 
any perturbation of a quantity can be expressed as the superposition of normal modes, we can decompose the spatial 
part into the summation over a series of harmonics Q*-''-' which are covariantly constant Q*-*^-* — and are chosen to 
be the eigenfunctions of the Laplace-Beltrami operator 

02g(fc) ^ _^Qik) . (30) 

We can then define a comoving scale A = 2Tra/k for each perturbation associated with a harmonic function Q^'^\ In 
our application, the harmonic decomposition is particularly useful as it allows us to distinguish the specific situation 
where the wavelengths of the perturbations are much larger than the Hubble scale (27ra/fc ^ H"^), in which case the 
Laplacian operator in equation (|30|l that is proportional to /c^, can be eliminated, yielding differential equations which 
are easier to solve. Although the use of a plane wave description is mathematically incorrect in curved space, it turns 
out that the only difference that arises is the allowable values of the wavenumbers. For a ffat geometry {K = 0) the 
eigenvalues form a continuous spectrum where A:^ > 0. The spectrum for an open model {K = +1) is discrete with 
k- = a{a + 2) where a = 1 . 2 , 3... A spacetime with negative curvature (K — —1) can accommodate the eigenvalues 
k^ = l + a'^ where > q 
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B. Governing equation for the interaction variable T°' 

In section (B.l) of Bctschart et al. 9], the solution for /" obtained by solving its wave equation was found to agree 
with the result calculated from the multiplication of the time dependencies of the shear and background magnetic 
field, determined individually. Solving the wave equation for the interaction variable, however, requires a harmonic 
decomposition. Since the interaction variable is not necessarily divergence-free, D°/a — cjab^'^B^ 7^ in general, 
it has a non-zero scalar contribution. This contribution is however exactly equal to Dq/?" and therefore drops out 
of Maxwell's equation Ijl^^ll rendering it purely solenoidal. Although decomposing the interaction term X° as a pure 
vector (as in the homogeneous case) is therefore incorrect, it is still possible to find the solution for X": either by 
solving for aab and the first-order part of B'^ separately and then multiplying the solutions obtained, or by proceeding 
in a similar manner as in 

GWs are purely tensorial and so we expand the representative shear variable with the help of tensor harmonics 

a., = ^aWQ«, (31) 

k 

where as usual Q^Qf,> = and D^Q^j^-* = ^{^'^ 1 0''^)Q^ab hold. Each gravitational wave mode is associated with the 
physical wavelength 

Agw = 27ra/fc . (32) 
The expansion of the magnetic field in pure vector (solenoidal) harmonics is 

B, = ^B(")g(") , (33) 

n 

and these vector harmonics obey the relations (9<a> = and D^Qi"' = —{n^/a'^)Q^J^\ Similarly as above, we 
associate with a given wavenumber n characterising a magnetic perturbation a characteristic length scale, 

= 27ra/n , (34) 

which we relate to the size of the magnetised region. 

To simplify the treatment of the interaction between GWs and the magnetic field, we proceed as follows. First, 
we assume the magnetised plasma region has a finite size = 27ra/m corresponding to some wavenumber m, which 
encodes the magnetic inhomogeneity over this region. We therefore write Ba = X]n>m>o 

bmce 

the interaction is most effective if the gravitational wavelength Aqw matches the size of the magnetic field region 
Ag 0, we may restrict ourselves to the resonant case where the gravitational and magnetic wavenumbers agree. 
This means that the comoving scale of the magnetic field perturbation is the same as that of the GWs, i.e. k = m. 
Consequently, the m-mode of the shear is the main contribution to the interaction, which now reduces to a single 
mode-mode term: la » ^("'^^('"'Qi^^Q^^). Making use of these considerations, it is now a straightforward task to 
obtain a closed equation for the interaction variable Xq by combining, as in , the standard evolution equations for 
the shear and the electric Weyl tensor, and using -B<a> — 0. In this way we readily arrive at 



i<a> + ^Xa + |eX<a> + (^O^ + - |p + |A) X„ - . (35) 
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Here, the second term on the LHS stems from the expression (— D^daf,) which emerges during the calculation. As 
expected, equation H35|) is equivalent to multiplying the m-mode shear equation (|12|) with the magnetic field B"^. We 
remind the reader that the derivation of equation (|35(l does not rely on the MHD approximation and is indeed valid 
for all values of the curvature index K and is also independent of the equation of state. 

C. Governing equation for T"" 

In the ideal MHD limit, the electric field is expressed as the cross-product of the primary magnetic field and the 
velocity, Ea = —eabcV^B'^- Using the linear velocity propagation equation together with the linear evolution 
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equation JSJ for the magnetic field, one easily finds E<a> + ©£"0 = 0, and therefore 

(curl£;a)^ = curl£;<Q> - |ecurl£;Q 

= -|ecurlK • (36) 

This further implies that the evolution equation of the MHD contribution encoded in !F°' is simply given by 

^<a> + f = . (37) 

It follows that the term Ta evolves like the electric field E'q, and decays as ~ a~^. 

VI. SOLUTIONS FOR SPATIALLY FLAT UNIVERSES 

Because the derived MHD equations above are only valid in the cold plasma limit, we only investigate the dust 
solutions for spatially flat models with zero cosmological constant. It is convenient to use the dimensionless time 
variable r = |i?o (^ — ^o) + 1 introduced in 0. In terms of this variable, the Hubble parameter evolves simply as 
Hq/t and the scale factor obeys a = aoT^/"^, where the zero index indicates evaluation at some arbitrary initial time 

to. 

Employing the time variable r, the equation for the interaction variable H35|l may be written as 



_ T " -A T ' 



/ TO 



2t2 V«oi^o 

while the equation H37|) for the MHD term transforms into 



r 3 



(38) 



whose solution is simply 



<a> I 3-7- ~ ' 



(39) 



(40) 



The generated magnetic field will typically depend upon x = m/{aoHo) = 27r(AH/A^)o = 27r(AH/AGw)o, which is the 
ratio between the size of the magnetised field region and the horizon size when the interaction begins. 

A. Limiting case where a; — > 

If the value of x is so small that we can drop the term x'^ in equation 1)38(1 . the general solution for the interaction 
variable is 



Ja{r)=C\T-s +C2T- 



(41) 



where Ci and C2 are integration constants. Since for the generation of magnetic fields the dominant mode is more 
important, we set C2 — and obtain the solution 



X,(r)=X°r~^=aS,Bgr-i. 
It is now a very simple exercise to integrate the induction equation H29() to find 



(42) 



(43) 



where the integration constant was determined by requiring the generated magnetic field to vanish initially. It follows 
that the total magnetic field measured by the fundamental observer due to the zero-zero mode interaction becomes 



(0) , 



Hq I Qq 



2 T 



(0) 



3 HoBi'^ 



(44) 
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where the second term in the square bracket originates from the interaction of the magnetic field Ba with the shear 
(already obtained in 0), while the third term represents the MHD contribution, which is the interaction of the 
magnetic field Ba with the plasma velocity perturbation Va- Note that the MHD contribution slowly decays away 
as the scale factor increases, in contrast to the gravito-magnetic part, which linearly grows with the scale factor. 
However, in the long-wavelength limit the contribution due to the GWs is negligible since one typically finds the 
shear anisotropy to be very small, (a/H)^ <^ 1. Moreover, since we may approximate J-o = Fq ~ (wi?/Ag)o due 
to Fa = —curlEa and the MHD relation H19|l . we see that the factor !Fq/{HqBq) in (|44|l is proportional to a; <C 1, 
and the MHD contribution turns also out to be negligible. Consequently, in the long-wavelength limit, there is no 
amplification of the initial magnetic field by GWs or velocity perturbations and one needs to consider the general case 
in order to look for an amplification. 



B. General case with a; 7^ 



When X is not negligible (the magnetised region is strictly finite), the general solution to the equation (|38|l is found 
to be 



Ia(r) 



DiJi 



(45) 



where Di, D2 are integration constants and Ji, J2 denote Bessel functions of the first and second kind, respectively. 
Since we are only interested in the dominant contribution, we set D2 = as before, noting that the Bessel function of 
the second kind is decaying on super-horizon scales a; ^ 1. The remaining integration constant takes then the value 



Di = 



4x^sin(2a;) — 3sin(2a;) -|- 6a;cos(2a;) 



Co So 



(46) 



Assuming that the induced magnetic field is zero initially when the interaction begins (t = 1), the solution for the 
rescaled magnetic field then becomes 



B{t) = 



erg Bp 
Hoy 



i sin(2a;) — xcos(2a;) -I- a; cos(2a;r 3 )r ^ — | sin(2a:r3 )r 



2£b 
3 Ho 



- 11^ T-i - 1 



(47) 



where we defined y = 4a;^ sin(2a;) — 3sin(2a;) + 6a;cos(2a;) and made use of (|39|l . Notice that in the limit a: — > one 
recovers the result Had we instead focused on the other branch of the solution we would find l|T7|) again 

but with the sin and cos functions as well as some signs interchanged. Hence, the m — m mode interaction generated 
magnetic field as seen by the fundamental observer moving with 4-velocity u° is given by the expression 



(m) 



i sin(2a;) 



a;cos(2a;) 



2 F 



(m) 



Ho ' 4a;2 sin(2a;) — 3sin(2a;) -|- 6a;cos(2a;) 3 HqBq 



(m) 



0(0" 



(48) 



here, the non-displayed terms are decaying with time and therefore irrelevant for the amplification process but can 
be inferred easily from H47|l if required. 

What happens if we take the full solution (|45|l of the interaction variable into account instead of just looking at 
one branch? Since both the Bessel functions of the first as well as the second kind are merely sin and cos functions 
modified by the same damping envelopes, one should in principle consider both branches even though their asymptotic 
behaviour is different. In this case, the integration constants Di and D2 are determined by requiring that initially 

_^^/(m)^(m) rpj^^ exact solution for the rescaled magnetic 



one has 2o = = 1) = ctq^^ Bo"^'' together with Xq' = a, 



field is then found in analogy with the above example [cf. H47|) ]. but it is too large to display it in full. However, the 
dominant contribution to the generated magnetic field can be written down in simple terms as 



B 



("-"'(a) ^B^'") (^) 



3 2a^)-Ha, 



2x2 



2 F 



(m) 







Hn 



3 HoB^o 



(m; 



(49) 



where it was again assumed that there is no generated magnetic field initially. If we employ the natural length scale 
Ag, inherent to the problem under consideration, we may estimate cto ~ (o'/Ag)o, implying CTq' « 2/3(<T/iJoAg)o, 
and also Fq « (u_B/Ag)o due to Fa = — curlE'a and the MHD relation H19() . Remembering further the definition of 
X = 2tt{Xyi/ Xg)o, we can finally write down the expression for the total magnetic field in a more convenient form: 



47r2 I Ah 



Ho 



Ah 

-^B/ 



Ah 

Afi 



0{a- 



(50) 
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This is our main result - it shows in detail how the magnetic field, resulting from the interaction of the background 
magnetic field Bq with GWs and velocity perturbations Vq in the plasma, depends on the initial conditions. 

Observe that at super-Hubble scales the MHD contribution becomes completely negligible mirroring the observation 
that plasma effects are typically more important on small scales. Our main result directly generalises our previous 
result (49) in 9J derived for the case of a homogeneous magnetic field to the inhomogcneous case. It should be stressed 
that the use of ideal MHD in the cold plasma limit allowed for a self-consistent treatment of the electric fields and 
plasma currents. 



VII. DISCUSSION 



If we look only at super-horizon scales and divide the result H5()|l through the energy density of the background 
radiation, /i^, (which decays in the same manner as the original magnetic field), the dominant contribution is then 
given by 



B 

~T/2 



a 



10 Vah;o ^^^0 




(51) 



where the wavenumber indices have been suppressed and the zero suffix indicates the time when the interaction begins. 
This result l|51|) was already found in Betschart et al. and previously reported by Tsagas et al. 10], a paper which 
employed the weak field approximation. The result H51|l can be applied to the reheating phase of the Universe at the 
end of inflation, for which the effective equation of state was that of dust (cf . |^ 0| for an application) . 
On the other hand, at sub-horizon scales the main part of the magnetic field is given by 



B 

~Tj2 



1^2 Ah 




(52) 



which could be applied to the matter-dominated phase of the Universe. In order to obtain an order-of-magnitude 
estimate we assume that the velocity perturbations (resulting from Thomson scattering) in the primordial plasma in 
effect start to interact with the pre-existing magnetic field B somewhat after matter-radiation equality with a redshift 
of roughly (zeq — 10^). The horizon at matter-radiation decoupling was = H^^{\ -\- Zoq)~'^/^ ~ 10~^ Mpc, where 
Hq denotes today's Hubble constant. A typical size of a seed field required for the dynamo mechanism is ~ 10 kpc on 
a comoving scale today, hence A^^** = (1 + ^cq)"^ A^ ~ 1 pc. From cosmic microwave background (CMB) measurements 
we know that at decoupling (zdoc — 10'^) the velocity perturbations satisfy Udcc < 10~^, while we read off from equation 
(|27|l that they decay like a^^. It follows that Vdec = i'eq(aeq/a) = v^qT^"^/^ = (1 + Zdoc)(l + 2;cq)^^fcq — 10~"'^t'eq, 
and whence Uoq — 10"'*. Taking everything together, one obtains the boost factor in (|52|l to be of order unity, that 
is, ('i'AH/Ag)^q ~ 1. Given that we can rely on the MHD approximation at the early stages of the matter-dominated 
scenario, the velocity perturbations in the plasma will lead at best to a doubling of the initial magnetic field strength. 

Comparing the result (|51|l with the final solution presented in Betschart et al. indicates that the magnitude of the 
amplification due to the interaction between GWs and magnetic fields is proportional to the square of the ratio of the 
coherence length Ag of the initial magnetic field and the initial size of the horizon Ah in both the homogeneous and 
inhomogeneous magnetic field cases. The additional MHD part of the field in H50|l above arises from the forcing term 
F'^ = — curl i?a = 2D'' [v[aBb]) whose time behaviour is obtained from a first-order propagation equation in which a 
Laplacian does not appear, its general solution therefore being independent of wavenumber and scale (of course, the 
initial conditions still depend on the size of the interaction region). It was found that the seed field's interaction with 
GWs is only important at super-horizon scales, while the interaction with plasma velocity perturbations dominates 
at sub-horizon scales. It is worth pointing out that there is no amplification at all in the long-wavelength limit. 

In contrast with the results in Betschart et al. , the generated magnetic field modes now have wavenumbers that 
are constructed from those of the interacting field and GWs and thus differ from those of the GWs alone. If the GW has 
wave vector ka and the 'background' magnetic field a wave vector (where Ag = li^ajm corresponds to the size of the 
magnetic inhomogeneity) , the wavenumber I of the induced field satisfies = {ka + ma){k°' + m°') = k^ + + 2mak'^. 
If for simplicity the first order wave vectors are assumed to be orthogonal, k°'ma — 0, then the m — m interaction 
yields a field with wavenumber £^ = 2m^ (4to^ in the parallel case). It follows that the characteristic wavelength 
A = 27ra/^ of the induced field is somewhat shorter than the original B-field by the superposition of the corresponding 
magnetic and GW wavenumbers. 

It is worth noting that the presence of a spatially homogeneous field is only consistent in the flat Universe 
whereas no restrictions on the spatial geometry arise at any stage in the derivation of the evolution equations with 
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(D''S°) ^ in this paper. This makes sense when we consider the spectrum of allowed wavenumbers for different 
geometries. For open models {K = +1), the lowest wavenumber is n = 3. For closed models {K — —1), we find that 
n > 1. However, the spectrum of wavenumbers in a flat spacetime {K — 0) is continuous, with n > 0. Given that 
the spatial homogeneity of the background field in restricts its associated wavenumber to n = 0, we see that this 
eigenvalue can only be accommodated in a flat Universe. 

More interesting is the relationship between curvature and magnetic fields. Einstein's th eory is geometrical which 
implies that vectors are directly coupled to the spacetime curvature via the Ricci identity |2g|. In [23, Tsagas and 
Maartens find that the evolution equation of the spatial gradient of the magnetic field contains a term CacdB^H^f ^ 
0{tB(-g)i which alludes to a non-local coupling to curvature. This result is confirmed in |^ with the appearance of a 
term containing the spatially projected Riemann curvature tensor in the propagation equation of (DaSf,) and indicates 
that the curvature sources magnetic inhomogencities. In the case of a homogeneous magnetic field in Bctschart et 
al., the spatial gradients of the magnetic field are at least second-order. In order to preserve the spatial uniformity 
of the first-order field through time, the spatially projected Riemann tensor may have to vanish to prevent it from 
sourcing (DaSb) so that the spatial gradients remain small and continue to contribute at higher-order only. In the 
analysis of the inhomogeneous magnetic field presented here, magnetic spatial eddies exist at first-order and for this 
reason, the boost from the coupling between the field and the curvature need not necessarily be eliminated. In fact, 
such couplings are explicitly retained in our approach via the standard commutation relations (cf. the appendix). 

VIII. CONCLUSION 

Although the focus of magnetogenesis in recent years has been the generation of large-scale magnetic fields, a self- 
sufficient mechanism still evades us. The galactic dynamo is indeed physically feasible and has been shown to generate 
fields with strengths matching current observations, but requires a reasonably strong seed field to work. To make 
this theory more robust, we need to find a way of producing seed fields that are suitable for subsequent amplification 
by the dynamo. Cosmological perturbations have been identified as a possible source of primordial magnetic field 
amplification which is present in the pre-recombination era. Considering second-order couplings between electrons, 
photons and protons, the electric current induced by the plasma vorticity (a known source of magnetic fields) and the 
additional contribution of the photon anisotropic stress are found to yield a magnetic field that is a sufRcient seed for 
the dynamo to work ^^S] ■ The coupling of density and velocity perturbations that are naturally occurrin g in the early 
Universe are shown to lead to similar resultant fields in the context of a relativistic charged multi-fluid |29l |. 

In this paper we built on previous work in Betschart et al, in which the full set of equations determining the 
evolution of the gravitational waves and the generated electromagnetic fields was presented, initially for the case of 
a homogeneous magnetic field, and generalised the analysis to the case of a spatially inhomogeneous magnetic field 
using the magnetohydrodynamic approximation (restricting ourselves to the dust case). Analysing the equations for 
a spatially flat dust FLRW Universe, we were able to confirm our previous results. In particular, we do not find any 
amplification in the case of the long- wavelength limit. 

Over and above the presentation of a physically viable mechanism for primordial magnetogenesis, this paper also 
establishes a formalism which guides the choice of proper second-order gauge-invariant variables. Using this method- 
ology, one is able to obtain results in terms of clearly defined quantities, with no ambiguity concerning the physical 
validity of the variables. 

The possibility of obtaining seed magnetic fields of sufRcient magnitude via a combination of inflationary physics 
and standard MHD theory is an exciting prospect. However, many of the parameters involved in estimating the size 
of the effect are not well known. For instance the spectrum of gravitational waves predicted by inflationary theory 
still has to be verified, as well as the large-scale structure of cosmological magnetic fields. For this purpose, more 
precise and complete measurements (such as the Planck mission) of the, e.g. polarization of the CMB would give 
much needed information js^f. Moreover, studies of interactions between incoherent gravitational wave distributions 
and turbulent magnetic fields could also be done in order to obtain a more detailed picture. This is left for future 
research. 
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APPENDIX: COMMUTATION RELATIONS 



Here we present various commutator relations which have been used in the text. The relations are given up to second 

order in our perturbation scheme. The vanishing of vorticity. LUab = 0, is assumed throughout in conjunction with the 
constraints Dq/x = D^p = which isolate the pure tensor modes. All appearing tensors are PSTF, Sab = S^ab>^ and 
all vectors Va, Wa are purely spatial. 
Commutators for first-order vectors Va' 

{T>aVb)^ = BaVb - iQBaVb - ^DeH + HJ Cdbc V (A.l) 

(curl Va)^ = curl K - curl K - aj' curl 14 - Hab (A.2) 
D[a ^b]Vc = [|e2 - i (m + A)] V^a hb]c + (^B sigma^[a - E^la) H] 

+hc[u{Eb]d-k^^b]d)V^ (A.3) 

Commutators for first-order tensors Sab- 

{BaSbc)^_ = BaSbc - lOBaSbc - a^BaSbc + 2HJ e^eib (A.4) 

[H^'Sab)^ = D^5„6 - |eD*5„6 - a'^B.Sab + eabc if'd S'"' (A.5) 

(curl^afc)'^ = curl 5a6-^e curl 5„6-(Te'ecd(a £''56)'* + 3^?c<a56> (A.6) 
curl curl 5a6 = -£>^5„6 + (/x + A - ^6^) 5„6 + i£'<a£>"56>c 

+35e<„ {E,^' - (A.7) 

Commutators for second-order vectors Wa'- 

{T)aWb)^ = T)aWb-\QT)aWb (A.8) 

D[„ T)b]Wc = [^e^ - I (M + A)] W^a hb]c (A.9) 
curl curl W„ = -D^W^ + D„ (dlvVF) + § (/x + A - ie') W„ (A.IO) 
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